Resonant Fermi gases with a large effective range 
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We calculate the equation of state of a Fermi gas with resonant interactions when the effective 
range is appreciable. Using an effective field theory for large scattering length and large effective 
range, we show how calculations in this regime become tractable. Our results are model-independent, 
and as an application, we study the neutron matter equation of state at low densities of astrophysical 
interest 0.002 fm -3 < p < 0.02 fm -3 , for which the interparticle separation is comparable to the 
effective range. We compare our simple results with those of conventional many-body calculations. 
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The properties of dilute Fermi gases when all two-body 
scattering lengths a s are large compared to the interpar- 
ticle separation r s and the range of the interaction R is 
small compared with r s are universal, in the sense that 
properties do not depend on details of the interparticle in- 
teraction Q . This is because the only dimensionful scale 
is the Fermi momentum fa, and the corresponding en- 
ergy scale is the Fermi energy £p = fcp/(2m), m being the 
fermion mass. Consequently all macroscopic observables 
are given by appropriate powers of fep or ep multiplied by 
universal factors. For instance, for two spin states with 
equal populations the energy per particle E/N of cold 
gases of, e.g., 6 Li, 40 K atoms or neutrons under these 
conditions is given by 
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where the universal factor £ is a number. 

Resonant, dilute Fermi gases were realized for the first 
time by O'Hara et al. in 2002 2]. The factor £ has been 
determined experimentally by a number of methods 0>S 
0, in atomic gases where the interaction may be tuned 
by controlling the magnetic field. The extraction of the 
equation of state leads to £ = 0.51 ± 0.04 0, £ ~ 0.7 0] 
and £ = 0.27 ^Q Qg 5], for temperatures in units of the 
Fermi temperature T/T F « 0.05 (except for T/T F » 0.6 
in 0). To date, the most reliable theoretical results are 
from T = 0, fixed-node Green's function Monte Carlo 
simulations, £ = 0.44 ± 0.01 6] and £ = 0.42 ± 0.01 0. 

The purpose of this letter is to consider the neutron 
gas as an example of a Fermi system with a large scat- 
tering length. For neutrons, the scattering length is 
unnaturally large and the currently accepted value is 
a nn = —18.5 ± 0.3 fm (for a recent review on the ex- 
perimental situation see ||). This is to be compared 
with the range of nuclear interactions, which is given 
by the mass of the lightest exchange particle, the pion, 
and thus R ~ l/m* ~ 1.4 fm. The effective range r c is 
expected to be approximately charge independent, and 
therefore we take the neutron-neutron effective range to 
be r nn = 2.7 fm, the same as the neutron-proton one 
Consequently, the neutron effective range is also signifi- 
cant, with TeTOTr ~ 2. 



For neutrons, the validity of the universal equation of 
state applies to densities with kpr c <C 1. This restricts 
neutron densities to p = kp/(3ir 2 ) < 10 -4 frn - , which is 
just below neutron drip density p n< j = 2.3 • 10 _4 fm -3 ~ 
10~ 3 /Oo at which neutrons become unbound. Here po — 
0.16 fm -3 is the saturation density of symmetric nuclear 
matter. Therefore, it is important to generalize the equa- 
tion of state to resonant Fermi gases with an appreciable 
effective range. In the regime kpr e ~ 1, the energy per 
particle can be expressed as in Eq. Q but with a system- 
specific factor £(k-pr c ), 
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This is the case for neutron matter at subnuclear den- 
sities p < po/lO. For the extension Eq. Q to be valid, 
the contributions to the energy from higher partial waves 
and many-body forces must be negligible. This ensures 
that neutron densities with kpr c ~ 1 do not probe further 
details of nuclear forces. In addition, the effective range 
expansion truncated at r e , kcotS(k) — —l/a s + r e k 2 /2, 
must describe the phase shifts 6(k) for S-wave scattering 
up to relative momenta k ~ k-p. The extension of the uni- 
versal energy to kpr c ~ 1 is also valid for trapped Fermi 
gases where both the scattering length and the effective 
range are tuned to large values. In this letter, we calcu- 
late £,(k-pr c ) for resonant Fermi gases with kpr c > 1 and 
explain how this regime becomes theoretically tractable. 

For neutron matter, the non-S-wave contributions 
to the energy per particle can be estimated from the 
Hartree-Fock results shown in Fig. The Hartree-Fock 
calculations are based on the model-independent two- 
nucleon interaction V\ ow k for different values of the mo- 
mentum cutoff A. The evolution of Mow k with the cutoff 
follows a renormalization group (RG) equation that guar- 
antees that nucleon-nucleon phase shifts are cutoff inde- 
pendent. We have previously found that the evolution of 
all large-cutoff nuclear interactions to A < 2.1 fm -1 leads 
to the same low-momentum interaction V\ ow & |l(| . Vary- 
ing the cutoff gives an approximate measure of higher- 
order contributions beyond Hartree-Fock and the effects 
of omitted many-body forces, since both are needed to 
obtain cutoff-independent results. We find in Fig.Q]that 
the low-density equation of state is cutoff-independent 
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FIG. 1: Hartree-Fock results for the energy per particle of 
neutron matter calculated from V\ ov , k for different cutoffs. 



for kp < 1.0 fin - , with negligible I 1 contributions 
for kp < 0.8 fm -1 . In addition, the effective range ex- 
pansion describes the 1 Sq phase shifts well up to relative 
momenta k < 1.0 fin - , with a deviation of —0.4, 0.2 and 
2.1 degrees from the empirical phase shifts 57.7, 48.6 and 
39.5 degrees for k — 0.4, 0.6 and 0.8 fm -1 respectively. 
Consequently, the equation of state of neutron matter is 
given by Eq. f° r ^f^c < 2. 

We calculate the equation of state Eq. J2J) using an 
effective field theory (EFT) for large scattering length 
and large effective range. EFT offers a systematic ap- 



proach to interactions at low energies, and therefore our 
results are model-independent. Under the conditions of 
interest, both the scattering length and effective range 
are low-momentum scales l/a s ~ l/r c ~ Q. The cor- 
responding EFT is realized by introducing a di-fermion 
field d |ll| (for a review see Ref. 01) with the lowest- 
order Lagrangian density given by (in units h = m = 1) 
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Here tp denotes the fermion field, and A and g are low- 
energy constants which describe the propagation of the 
di-fermion field and its coupling to two fermions respec- 
tively. We also note that ido + V 2 /4 is the operator 
for the two-body energy corrected for the center-of-mass 
motion. The EFT Lagrangian Eq. @ was introduced for 
nucleon-nucleon scattering in Ref. 11] following work due 
to Weinberg In the di-fermion EFT, the T matrix 
for two-body scattering depends only on the energy E of 
the interacting particles in the center-of-mass system, 



T(E) 



A-E-g 2 I (E + ie) 
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with J (.E) = Jq d 3 p/{2n) 3 (E-p 2 )- 1 . This corresponds 
to summing diagrams where the di-fermion propagator 
(represented as a double line) is dressed by fermion loops, 
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Iq(E) denotes the loop integral regulated by a momen- 
tum cutoff A and the on-shcll T matrix is given by 
4ttT^ 1 (E = k 2 ) = l/a s — r c k 2 /2 + ik, which can equiv- 
alently be obtained from an energy-dependent potential, 
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which follows after integrating out the di-fermion field. 
Matching the lowest-order low-energy constants A and g 
to the effective range expansion yields 
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Consequently, for systems with large scattering length 
and large effective range the potential V(E) scales at low 
energies E - Q 2 as V{E) - (Q+A+Q 2 /A)-\ The cutoff 
generates higher-order terms ~ fc 4 /A 3 in T^ 1 (E = k 2 ), 



which are suppressed for fc/A <C 1. Therefore, we con- 
sider the lowest-order EFT with large cutoffs. This is 
not required for V| ow k, since the RG generates all higher- 
order contact interactions necessary to maintain cutoff- 
independent two-body observables. If one rescales the 
di-fermion field by d — > g d, the interaction terms be- 
come coupling-independent and the kinetic term acquires 
a factor 1/g 2 . For large cutoffs 1/g 2 ~ r c , and therefore 
the di-fermion field is an auxiliary field for positive effec- 
tive range, even if there is a two-body bound state with 
a s > 0. Integrating out the d field gives the standard 
contact interaction EFT. For negative effective range the 
di-fermion kinetic term has the normal sign and the La- 
grangian Eq. © is the atom-molecule model. The great 
advantage of the present formulation is that it enables 
one to study systems with either positive or negative r e . 

Next, we study the scaling of diagrammatic contribu- 
tions to the energy. Particle-particle loops are restricted 
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only by the cutoff, and therefore these scale as ~ A 
for large cutoffs. By contrast, particle-hole and hole- 
hole intermediate states scale with the low-momentum 
scale ~ Q. Particle-hole loops, for instance, scale with 
the density of states ~ kp. Since the lowest-order inter- 
action scales as V{E) ~ (Q + A + Q 2 /A) -1 , cutoff inde- 



pendence for large cutoffs minimally requires summing 
particle-particle ladders with Pauli blocking, such that 
the A -1 scaling of the potential cancels with each loop 
~ A 0| . Neglecting terms which vanish for large cutoffs, 
we find for T me d, the on-shell in- medium T matrix 
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where P denotes the particle-pair momentum and we 
have neglected the imaginary part, as the latter van- 
ishes for k < y/k F — P 2 /4 and therefore does not con- 
tribute to the energy. The resulting in-medium T ma- 
trix can be viewed as an effective interaction and scales 
as T me d(E — k 2 ;P) ~ Q^ 1 independent of the cutoff. 
In addition to the scattering length and effective range 
terms, the inverse effective interaction acquires a density- 
dependent contribution — (kp + P/2)/w and two loga- 
rithmic terms. The first logarithmic term includes the 
BCS singularity for P = and k — k-p. The result- 
ing contribution to the energy is integrable and gener- 
ally small. For neutron matter, realistic calculations of 
the superfluid pairing gap including polarization effects 
give A/ep ~ 0.1 for densities with kpr ~ 1 We 
therefore expect corrections to the energy per particle 
£si{kpr c ) ps —5/8 (A/ep) 2 ~ —0.01 due to pairing in this 
regime. The second logarithmic term in Eq. JBJ is regu- 
lar, and for P = it simplifies to —kp/n. 

For large effective range, interactions are weaker at 
higher momenta and the effect on the in-medium T ma- 
trix can be understood qualitatively by inserting average 
values for the momenta k 2 = 3/10 fc F and P 2 = 6/5 k F 
into Eq. JHJ. With this qualitative estimate, the average 
effective interaction for resonant Fermi gases is given by 
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with Uq = 0.56. In the universal regime with kpr D <C 
1, the average effective interaction is given by T me d = 
— Air/(Uokp). The remaining contributions to the energy 
are of particle-hole or hole-hole nature, and scale as Q ~ 
kp. In this EFT it is thus necessary to sum all diagrams, 
as the effects of each additional vertex ~ fcp 1 combined 
with the additional loop ~ kp are of order one. 

For Fermi gases with large effective range, the loop 
scaling is unchanged but additional T mec j vertices are now 
suppressed by (1+0.27 /cp^c) -1 - This implies that correc- 
tions to the Hartree-Fock equation of state with T me d as 
effective interaction come with powers of (1+0.27 fc^e) -1 
relative to the universal case. In perturbation theory, 
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FIG. 2: Result for ^hf(^f^c) (solid line) which describes the 
equation of state of resonant Fermi gases with large effective 
range fere > 1. The dashed lines provide an error estimate 
for £hf as discussed in the text. 



particle-particle-irreducible contributions to the energy 
enter in third order. With T mc d = 4-7T a me( j , we esti- 
mate the leading non-ladder corrections to the £(kpr e ) 
obtained from T me( j at the Hartree-Fock level as A£ = 
0.077 (U (l + 0.27fc F r o ))- 3 _16]. This only provides a 
qualitative estimate, since we evaluated the in-medium 
T matrix for average particle configurations. 

To a good approximation, the equation of state of res- 
onant Fermi gases with kpr c > 1 can therefore be calcu- 
lated from T me d at the Hartree-Fock level, 
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with l/a s — in Eq. JSJ. Based on the arguments given 
above, we expect corrections to £ beyond Hartree-Fock 
of order A£ = 0.22 — 0.03 and £ s f ps —0.01 due to pairing 
for densities with kpr c = 1 — 5. Our results for ^hf^f^c) 
are shown in Fig. |3 with the theoretical error estimate, 
ft is intriguing that we find £hf(1) ~ 0.47, which is close 
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FIG. 3: Comparison of the energy per particle for neutron 
matter calculated from Eq. il I Uti for a s = —18.5 fm and r e = 
2.7 fm (solid line) to Fermi hypernetted-chain (triangles) (Tt| 
and Brueckner results (circles) (HI . We also show the V\ ov ,k 
Hartree-Fock results of Fig. The dashed lines provide an 
error estimate for £hf- 

to the universal factor £ = 0.44 ± 0.01 |6|. This implies a 
cancellation of higher-order contributions for kpr c <C 1. 

For kpr c ^> 1, the effective interaction is of the or- 
der T me( j ~ — (r c /cp) _1 and interaction effects disap- 
pear for infinite effective range. In the latter regime, 
£ w 1 + C/(k-pr c ), with a coefficient C < of order 
one. Consequently, £ — > 1 only for very large effective 
ranges. We find £hf(10) = 0.79 and £hf(100) = 0.96. If 
the effective range is negative, the average interaction in 
the medium is greater than for r c = 0, and higher-order 
terms are only reduced for very large \kpr c \ 3> 1. 

Finally, we take into account a large, but finite scatter- 
ing length, and present in Fig. [3] results for the equation 
of state of neutron matter in the regime kpr D ~ 1- This 
is straightforward by including the l/a s term in Eq. JSJ. 
Our results are model-independent and constrain neu- 
tron matter at subnuclear densities 0.002 fm -3 < p < 
0.02 fm~ 3 , which is of relevance to the physics of neu- 
tron stars and supernovae. We also compare our neu- 
tron matter results to many-body calculations using 
Fermi hypernetted-chain techniques or Brueckner 
theory 18] . As can be seen from Fig. all microscopic 
results lie within our theoretical error estimates. 

In summary, we have calculated the equation of state of 
resonant Fermi gases with large effective range kpr c ^ 1 
in the di-fermion EFT, with particular attention to neu- 
tron matter at subnuclear densities. In the lowest-order 
EFT, cutoff independence minimally requires the resum- 
mation of particle-particle ladders [13], which leads to an 
effective interaction that becomes weak for all kpr c ^ 1- 
The kpr c dependence was used to show how resonant 
Fermi gases with large effective range, such as low-density 
neutron matter, are theoretically tractable and to es- 
timate the error of our results. Neutron matter equa- 
tions of state obtained using conventional many-body ap- 
proaches are consistent with our model-independent di- 



fermion EFT results within these errors. For kpr c ~ 1, 
the energy is found to be close to the universal regime. 
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